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Characterisation of Culex quinquefasciatus (Diptera: Culicidae) 
larval habitats at ground level and temporal fluctuations 
of larval abundance in Cordoba, Argentina 

Marta GrechA, Paolo Sartor, Elizabet Estallo, Francisco Luduefia-Almeida, Walter Almiron 

Centre de Investigaciones Entomologicas de Cordoba, Institute de Investigaciones Biologicas y Tecnologicas, 
Consejo Nacional de Investigaciones Cientificas y Tecnicas, Universidad Nacional de Cordoba, Cordoba, Argentina 

The aims of this study were to characterise the ground-level larval habitats of the mosquito Culex quinquefascia- 
tus, to determine the relationships between habitat characteristics and larval abundance and to examine seasonal 
larval-stage variations in Cordoba city. Every two weeks for two years, 15 larval habitats (natural and artificial wa- 
ter bodies, including shallow wells, drains, retention ponds, canals and ditches) were visited and sampled for larval 
mosquitoes. Data regarding the water depth, temperature and pH, permanence, the presence of aquatic vegetation 
and the density of collected mosquito larvae were recorded. Data on the average air temperatures and accumulated 
precipitation during the 15 days prior to each sampling date were also obtained. Cx. quinquefasciatus larvae were 
collected throughout the study period and were generally most abundant in the summer season. Generalised linear 
mixed models indicated the average air temperature and presence of dicotyledonous aquatic vegetation as variables 
that served as important predictors of larval densities. Additionally, permanent breeding sites supported high larval 
densities. In Cordoba city and possibly in other highly populated cities at the same latitude with the same environ- 
mental conditions, control programs should focus on permanent larval habitats with aquatic vegetation during the 
early spring, when the Cx. quinquefasciatus population begins to increase. 
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Mosquitoes affect humans and animals by transmit- 
ting the disease-causing agents of several serious dis- 
eases. Such arboviral diseases likely represent the major 
threat to public health in Argentina, where the St. Louis 
encephalitis (SLEV) flavivirus has been isolated from 
Culex quinquefasciatus mosquitoes collected in central 
regions (Diaz et al. 2006). In 2005, the largest SLEV 
epidemic in South America and Argentina occurred in 
Cordoba Province, resulting in at least 47 human cases 
and nine deaths (Almiron et al. 2006, Spinsanti et al. 
2008). In 2010, 43 cases of SLEV were identified in Ar- 
gentina, four of which were reported in Cordoba Prov- 
ince (OPS 2010). 

The successful application of integrated mosquito 
control techniques requires an adequate understand- 
ing of mosquito larval habitats as well as the features 
of temporal dynamics. The number of studies address- 
ing natural and artificial mosquito larval habitats has 
increased in recent years in the central (Almiron & 
Brewer 1996, Fischer et al. 2000, Fischer & Schweig- 
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mann 2004, 2010, Pires & Gleiser 2010) and northeast- 
ern regions of Argentina (Stein et al. 2002a, b, 2005, 
2011). Cx. quinquefasciatus can use containers as larval 
habitats and shows the ecological plasticity to breed in 
both polluted and clear waters (Forattini 2002). Almiron 
and Brewer (1996) described and classified immature- 
stage mosquito habitats in Cordoba Province, where Cx. 
quinquefasciatus was grouped together with other Culex 
species that shared the same habitat requirements, with 
Cx. quinquefasciatus being found to be more frequent in 
artificial larval habitats. In another study, artificial and 
natural larval habitats (pools, ditches, canals, lagoons 
and water sources) were surveyed in Cordoba city dur- 
ing the summer months. Several species oiAedes, Culex 
and Mansonia were identified and Cx. quinquefasciatus 
was shown to be the most abundant and most widely dis- 
tributed in the city (Pires & Gleiser 2010). Furthermore, 
the temporal dynamics of Cx. quinquefasciatus popula- 
tions have been investigated in Buenos Aires (Fischer 
et al. 2000, Fischer & Schweigmann 2010) and Cordoba 
(Almiron & Brewer 1995), where this mosquito species 
was reported to be present throughout the year. 

Knowledge of mosquito larval habitats is essential 
for identifying the most productive sites in a given area 
and developing effective control methods. This study at- 
tempted to determine which ecological factors related 
to larval habitat dynamics drive Cx. quinquefasciatus 
abundance. The aim of this study was to characterise the 
ground-level mosquito larval habitats of Cx. quinquefas- 
ciatus, to determine the relationships between habitat 
characteristics and larval abundance and to study sea- 
sonal larval- stage variations in Cordoba city. 
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MATERIALS AND METHODS 

Study area - The study was conducted in Cordoba 
city (3r22'S 64°12'W), which is located at 360-480 m 
above sea level. Cordoba city is the main city and capital 
of Cordoba Province (Argentina) and has a resident hu- 
man population of 1,330,023 inhabitants (INDEC 2010). 
Cordoba city covers a land area of 576 Km^ and approxi- 
mately 37.2% of the city surface is urbanised, which is 
surrounded by agricultural fields (Jarsiin et al. 2003). 
The climate is temperate with distinctively hotter-wet 
(October-April mean maximum and minimum tempera- 
tures: 27.15°C and 14.87°C, respectively; mean monthly 
rainfall: 111.43 mm) and cooler-dry (May-September 
mean maximum and minimum temperatures: 19.66° and 
6.1°C, respectively; mean monthly rainfall: 17.98 mm) 
seasons (SMN 2011). 

Mosquito sampling - Cordoba city was partitioned 
into 100 plots of 4 km^ each, 52 of which were located in 
the urbanised area. These 52 plots were numbered and 
represented potential sampling positions. Fifteen plots 
were randomly selected in the urban public area. In each 
selected plot, preliminary samplings were performed to 
detect the mosquito larval habitats of Cx. quinquefascia- 
tus at ground level (natural and artificial water bodies, 
including shallow wells, drains, retention ponds, canals 
and ditches). Sampling plots near downtown were dis- 
carded because of the difficulty of finding these types of 
mosquito larval habitats and were replaced successively 
by other nearby plots to set up the scheme shown in Fig. 
1. Once detected, only one larval habitat within each 
plot was selected. Field sites were sampled for larval 
mosquitoes every two weeks from June 2008-June 2010 
(Fig. 1). Three samples were collected from each larval 
habitat using a standard white 350-mL dipper. The col- 
lected material was transferred to the laboratory in small 
plastic flasks. The third and fourth-instar larvae were 
killed and stored in 80% ethanol for taxonomic deter- 
mination according to available keys (Darsie 1985). The 
younger (1st and 2nd) instars were reared incubated until 
reaching the fourth instar in plastic trays with 750 mL of 
water from a natural larval habitat and were fed 0.25 mg 
of liver powder per larva daily. The water surface was 
cleaned each day with filter paper to avoid complications 
due to fungal and bacterial development. 

For each larval habitat and sampling date, the follow- 
ing characteristics were recorded: water depth (WD), 
temperature and pH, permanence, the presence of aquat- 
ic vegetation and the number of mosquito larvae collect- 
ed. Density was expressed as the number of larvae per 
dip (total number of larvae/number of dips). The aver- 
age WD was estimated from three measurements per- 
formed at random along the edges and in the middle of 
the larval habitat. Water temperature (WT) and pH were 
determined in the field with a digital thermo-pH-meter. 
Because the larval habitats were systematically sampled 
(the route, which included the selected larval habitats, 
began at the laboratory starting point at 08:00 am, ran 
through sites 1-15 and reached the laboratory end point 
at approximately 12:30 am, covering approximately 80 
km), the WT measurements were corrected to compare 
the obtained datasets using the following formula: 



rc., = ro, + ro.,-rp., 

where TC. ^ is the corrected WT from larval habitat /, reg- 
istered on date t, TO^ is the average WT of the 15 larval 
habitats on date t, TO. ^ is the WT from larval habitat /, 
measured in the field on date t, and TP. ^ is the predicted 
WT for larval habitat / on date t, calculated for each sam- 
pling date via linear regression between TO vs. the time 
of day (h) at which it was registered. Then, TP.Qi) = a + 
bh where TP.Qi) is the predicted WT from larval habi- 
tat ; and h is the time of day. The presence or absence 
(1/0) of aquatic emergent vegetation was recorded and 
the plants were broadly grouped into monocotyledons 
and dicotyledons. Water permanence was recorded as 
the percentage of sampling dates on which each larval 
habitat contained water. The daily air temperature and 
precipitation data for Cordoba city were provided by the 
National Meteorological Service. The average mean dai- 
ly air temperature and accumulated precipitation for the 
15 days prior to each sampling date were calculated. 

Data analysis - Data from the examined two-year pe- 
riod were analysed with a paired t test to compare the 
larval abundances pooled by date (paired between years). 
To evaluate the relationship between the observed Cx. 
quinquefasciatus larval density, habitat characteristics 
and meteorological variables, we used generalised linear 
mixed models (GLMMs) with a negative binomial family 
distribution and logarithmic link function (Crawley 2007, 
Zuur et al. 2009). The application of GLMM techniques 
allows ecologists to address elements such as nested data, 
temporal and spatial correlations, data with many zeros 
and heterogeneity of variance. Furthermore, negative 
binomial GLMMs are useful for modelling ecological 
datasets that show overdispersion, in which the variance 
is larger than the mean (Zuur et al. 2009). Such models 
are referred to as mixed because they include both fixed 
and random components. The fixed components are de- 
fined as the factors for which the levels are experimen- 
tally determined or the interest lies in the specific effects 
of each level, such as the effects of covariates or differ- 
ences among treatments and interactions, while the ran- 
dom components are defined as the factors for which the 
levels are sampled from a larger population or the interest 
lies in the variation among them rather than the specific 
effects of each level. In the present study, we applied a 
random intercept model wherein the intercept was al- 
lowed to change for each site. The WT, temperature and 
pH, presence of aquatic vegetation (monocotyledons and 
dicotyledons), average air temperature and accumulated 
precipitation 15 days prior to each sampling date were in- 
cluded as fixed effects. The larval habitat was included as 
a random effect. We chose the larval habitat as a random 
effect of the mixed model to account for the lack of in- 
dependence or pseudoreplication associated with the col- 
lection of multiple measurements and samples from the 
same sites over time. Additionally, we were not interest- 
ed in determining the exact nature of the site effect. We 
employed a backward selection procedure in which we 
removed non-significant terms from the model, one by 
one, in decreasing order of their p-values (Crawley 2007). 
Statistical analyses were performed with R software, ver- 
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sion 2.15.1 (glmmADMB package) (R Core Team 2012). 
To apply the mixed effects model (MEM) in R, we used 
the following code: 

> library(glmmADMB) > MEM<glmmadmb(datos$DL 
~T15+P15+pH+WT+WD+VAm+VAd+(l|Site), zeroln- 
flation = FALSE,data = datos,family = "nbinom") 
> summary(MEM) 

where the MEM was applied using the glmmadmb 
function. The response variable DL represents the lar- 
val density of Cx. quinquefasciatus. The included fixed 
effects were the averaged air temperature (TJ^) and ac- 
cumulated precipitation {P^^ on the 15 days prior to the 
sampling date as well as the pH, WT, WD and mono- 
cotyledonous and dicotyledonous aquatic emergent veg- 
etation (VAm, VAd). The random term IjSite specifies a 
random intercept model. The argument on the right hand 
side of the | symbol is a nominal variable (site). The sum- 
mary command provides the output of the model. 

RESULTS 

Characterisation of larval habitats - Six species be- 
longing to two genera were identified in the mosquito 
larval habitats. Cx. quinquefasciatus represented 99.64% 
of the total, while Culex hidens-mollis, Culex apicinus, 
Culex eduardoi-dolosus, Culex saltanensis and Ochle- 
rotatus albifasciatus represented 0.17%, 0.16%, 0.007%o, 
0.005%) and 0.02%) of the collected species, respectively. 
The average number of Cx. quinquefasciatus larvae ob- 
tained per dip varied from 0.04-263.6. The mean water 



pH ranged from 5.59-6.46, the WT from 14.95-22.2°C, 
the WD from 7.08-28.92 cm and permanence from 41- 
100% (Table I). 

Temporal variations in Cx. quinquefasciatus larval 
abundance - Cx. quinquefasciatus larvae were collected 
throughout the entire study period. The total number of 
larvae collected during the first season (59,410 larvae) 
was significantly greater (t = 2.66; p = 0.0151) than that in 
the second season (20.493 larvae). The initially recorded 
presence of immature stages gradually decreased during 
the winter period (June-September) in both seasons. How- 
ever, in early September 2008-2009 and October 2009- 
2010, the density of Cx quinquefasciatus larvae increased 
markedly. The highest larval density peak was recorded in 
February 2009 and another, less important, peak was ob- 
served at the end of March 2009. No considerable increas- 
es were detected during the second season (Fig. 2B). 

The mean air temperature and accumulated precipi- 
tation recorded during the first season were 19.53°C and 
563.4 mm, respectively. During the second study period, 
the mean temperature was 18.61°C and the accumulated 
precipitation was 1,085.2 mm. The average air tempera- 
tures and accumulated precipitation observed in the 15 
days prior to the sampling dates for the entire study pe- 
riod are included in Fig. 2A. 

Habitat characteristics and meteorological variables 
in relation to Cx. quinquefasciatus larval abundance - 
The variables from the GLMM that best described the 
relationship between the Cx. quinquefasciatus larval 




Fig. 1: map of Argentina. Location of 15 sampled larval habitats in Cordoba city (Cordoba Province), Argentina. 



Cx. quinquefasciatus larval habitats • Marta Grech et al. 775 



density and the examined habitat characteristics and me- 
teorological variables were the average air temperature 
during the 15 days prior to the sampling date (z = 2.54; 
p = 0.011) and the presence of dicotyledonous aquatic 
emergent vegetation (z = 2.47; p = 0.014). The other ex- 
planatory variables were not significantly related to the 
larval density (p > 0.05) (Table II). 

DISCUSSION 

In the present study, the average air temperature two 
weeks prior to each sampling date was found to be an 
important predictor of the Cx. quinquefasciatus larval 
density. However, the accumulated precipitation was not 
related to this variable, which could be due to the fact 
that almost all of the studied larval habitats with high 
larval densities were permanent and their water sources 
were mainly sewage and irrigation water, rather than 
rainfall. Relationships between larval abundances and 
meteorological variables have also been reported by oth- 
er authors, though with differing results being obtained. 
Almiron and Brewer (1995) observed that in artificial 
larval habitats in Cordoba city, which were mainly rain- 
filled, the larval abundance of Cx. quinquefasciatus was 
not correlated with variables such as air temperature and 
rainfall. However, in Buenos Aires city, the larval abun- 
dance of Culex pipiens was found to be positively associ- 
ated with the mean weekly air temperature (Fischer et al. 
2000, Fischer & Schweigmann 2004). 

The presence of dicotyledonous aquatic emergent veg- 
etation was also an important predictor of larval density 
in this study. Aquatic vegetation provides protection from 
predators and physical disturbances such as wind action 
and increases the availability of food, thus offering a more 



suitable habitat for immature mosquitoes (Russell 1999). 
Relationships between vegetation and larval abundance 
have been previously reported in Cordoba Province, 
where 42 localities were sampled to describe and classify 
mosquito larval habitats and the presence of aquatic veg- 
etation in Cx. quinquefasciatus larval habitats was indi- 
cated as the most important feature related to the larval 
abundance of this species (Almiron & Brewer 1996). Pires 
and Gleiser (2010) observed that, in natural and artificial 
larval habitats in Cordoba city, medium and high larval 
densities occurred in ditches and canals and were asso- 
ciated with the presence of vegetation. Positive associa- 
tions between increasing vegetation cover and immature 
Cx. quinquefasciatus abundance have also been reported 
in the Buenos Aires (Fischer & Schweigmann 2004) and 
Chaco Provinces (Stein et al. 2002a, 2005, 2011). 

Mosquito populations are sensitively dependent on 
the WTs in their larval habitats, which can directly affect 
larval survival and development times and indirectly in- 
crease decomposition processes and food availability 
(Clements 1992, Lindblade et al. 2000). In the present 
study, WT could not explain the observed variation in 
Cx. quinquefasciatus larval densities. However, the 
lack of an association between these variables does not 
represent sufficient evidence to conclude that the WT 
is less important in explaining larval density. The mea- 
surements of this variable conducted during our research 
only reflect the WT at a single point for each larval habi- 
tat and sampling date, rather than variations that would 
influence larval abundances. 

The mean pH values observed in this study were al- 
ways within the acidic pH range (range: 5.59-6.46) and 
were not related to the Cx. quinquefasciatus larval den- 



TABLE I 



Mean larval density (DL) ± standard deviation (SD) of Culex quinquefasciatus collected 
at 15 larval habitats and mean values of pH, water temperature (WT) and water depth (WD) registered in each site 





Type of 




Permanence' 


Aquatic 




WT 


WD 


Site" 


larval habitat 


DL±SD 


(%) 


vegetation 


pH 


(°C) 


(cm) 


2 


DR 


263.6 ±346.12 


100 


VAm, VAd 


6.01 ±0.39 


20.97 ± 4.42 


23.68 ± 7.43 


9 


CH 


131.48 ±240.47 


80 




6.00 ±0.38 


19.07 ± 4.99 


18.67 ± 10.49 


4 


CH 


123.66 ± 172.38 


100 


VAm, VAd 


5.89 ±0.54 


18.92 ±5.17 


19.15 ± 7.4 


13 


DR 


94.83 ± 118.4 


97 


VAm, VAd 


6.03 ±0.28 


18.46 ±5.6 


28.92 ± 11.87 


1 


SW 


74.29 ± 108.13 


100 




6.38 ±0.67 


17.36 ± 9.66 


9.66 ±2.95 


7 


CH 


73.46 ± 139.07 


100 


VAm, VAd 


6.02 ± 0.41 


18.57 ±4.75 


19.96 ± 5.24 


5 


CH 


58.14 ± 112.5 


100 


VAm 


6.01 ±0.5 


18.22 ± 5.42 


16.73 ± 6.25 


12 


CH 


21.67 ±39.79 


76 




5.96 ±0.12 


22.2 ±4.53 


11.21 ±3.96 


15 


SW 


14.18 ±24.37 


46 




6.44 ±0.51 


20.15 ±3.12 


7.31 ± 1.65 


11 


CH 


9.58 ± 29.23 


73 




6.00 ±0.34 


19.24 ± 5.62 


18.85 ± 13.4 


6 


CH 


2.98 ± 5.92 


65 




5.82 ±0.35 


19.39 ±4.24 


13.55 ±6.19 


14 


CH 


1.09 ±3.02 


73 


VAm, VAd 


6.05 ± 0.49 


19.53 ± 4.47 


12.66 ±6.98 


3 


RP 


0.74 ±1.75 


46 


VAm 


6.23 ±0.27 


17.83 ± 5.58 


8.39 ±3.59 


8 


SW 


0.11 ±0.27 


41 




6.46 ±0.4 


14.95 ± 4.69 


7.08 ± 3.83 


10 


CH 


0.04 ±0.15 


62 




5.59 ±0.64 


17.88 ±5.19 


17.24 ± 7.43 



a: sites are listed in decreasing order of mosquito density; b: percentage of sampling dates that each larval habitat contained 
water; CH: canals; DR: drains; RP: retention ponds; SW: shallow wells; VAm VAd: presence of monocotyledonous or dicotyle- 
donous aquatic emergent vegetation. 
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sity. Similar results have been obtained in other locali- 
ties in Cordoba Province, where the observed pH ranges 
were wider (from 6.4-8) and no relationship between this 
variable and Cx. qiiinquefasciatiis larval abundance was 
observed (Almiron & Brewer 1996). However, in Chaco 
Province, the presence of immature Cx. qiiinquefas- 
ciatiis stages was positively associated with pH values 
ranging from 6.8-8.5 (Stein et al. 2011). 

The presence of immature mosquito stages has been 
recorded at a wide range of breeding sites, from small and 
highly ephemeral to the large and permanent water bod- 
ies (Forattini 2002). Pires and Gleiser (2010) observed 
that the presence of Cx. quinqiiefasciatus larvae in Cor- 
doba city was related to temporary and semi-permanent 
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Fig. 2A, B: mean values of meteorological variables and mosquito den- 
sity during the study period (June 2008-May 2010). DL: larval density. 



TABLE II 

Variables from generalised linear mixed models 
explaining variation in Culex quinquefasciatus larval density 
with p < 0.05 and parameter estimates ± standard error (SE) 



Variable 


Parameter estimate ± SE 


z 


P 


T 

IS 


0.16 ±0.06 


2.54 


0.011 


VAd 


1.4 ±0.58 


2.47 


0.014 



T|j: average air temperature 15 days prior to each sampling 
date; VAd: dicotyledonous aquatic emergent vegetation. 



water bodies, both of which were associated with rain- 
fall as the water source. However, in this study, breed- 
ing sites with high water permanence values (80-100%) 
were found to support high Cx. quinquefasciatus larval 
densities. Because the eggs of this mosquito cannot tol- 
erate desiccation and hatching occurs when embryonic 
development is complete (Forattini 2002), permanent 
breeding sites might be more appropriate than tempo- 
rary water bodies for immature-stage development. 

In this study, the presence of Cx. quinquefasciatus 
larvae observed during the entire study period, even 
during the winter months, with a peak larval density be- 
ing recorded in February, suggests that the extreme tem- 
peratures that occur in Cordoba city (temperate climate) 
do not prevent larval development. The results of the 
present work correspond to other surveys of this species 
in Cordoba city, where Cx. quinquefasciatus adults and 
immature stages were collected during all seasons of the 
year and showed two peaks of larval abundance, in the 
spring and summer months (Almiron & Brewer 1995). In 
Buenos Aires Province (temperate climate), the presence 
of Cx. quinquefasciatus has also been reported through- 
out the year (Fischer et al. 2000, Fischer & Schweigmann 
2010) and immature stages of this species are generally 
most abundant in the summer (Fischer & Schweigmann 
2004). In subtropical Chaco Province, Cx. quinquefas- 
ciatus, together with Aedes aegypti, has been reported 
as the most abundant species in artificial larval habitats 
during the rainy and warm seasons; however, during dry 
and cold months, the Cx. quinquefasciatus larval abun- 
dance decreases in a manner consistent with the patterns 
observed in other temperate areas (Stein et al. 2002b). In 
the present study, the decreased larval density observed 
during the summer months of the second period could 
be a consequence of environmental sanitation measures 
(vegetation removal and larvicide application) applied 
during January and February 2009 to reduce mosquito 
larval habitats with the support of the local authorities 
(no available data). Other factors that were not exam- 
ined in this study, such as the presence of competitors or 
predators, might also affect larval abundance. 

More recent approaches to disease prevention and vec- 
tor control programs have focused on the cost-effective 
utilisation of limited resources to reduce the vector popu- 
lations (OPS 1995, Knight et al. 2003). Thus, in Cordoba 
city and perhaps in other localities at the same latitude 
and with the same environmental conditions, mosquito 
control programs should focus on permanent larval habi- 
tats with aquatic vegetation during the early spring, when 
the Cx. quinquefasciatus population begins to increase. 
Furthermore, the average air temperature during the 15 
days prior to collection could be used to develop predic- 
tive abundance models for this mosquito species. 
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